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Abstract Strong lensing is a powerful tool to address three major astrophysical issues: un- 
derstanding the spatial distribution of mass at kpc and sub-kpc scale, where baryons and dark 
matter interact to shape galaxies as we see them; determining the overall geometry, content, and 
kinematics of the universe; studying distant galaxies, black holes, and active nuclei that are too 
small or too faint to be resolved or detected with current instrumentation. After summarizing 
strong gravitational lensing fundamentals, I present a selection of recent important results. I 
conclude by discussing the exciting prospects of strong gravitational lensing in the next decade. 
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1 Introduction 

As photons from distant sources travel across the universe to reach our telescopes 
and detectors, their trajectories are perturbed by the inhomogeneous distribu- 
tion of matter. Most sources appear to us slightly displaced and distorted in 
comparison with the way they would appear in a perfectly homogeneous and 
isotropic unive rse. This phenomenon is called weak gravitational lensing (e.g. 



Refregier 20031 . and references therein). Under rare circumstances, the deflection 
caused by foreground mass overdensities such as galaxies, groups, and clusters 
is sufficiently large to create multiple images of the distant light source. This 
phenomenon is called strong gravitational lensing. Due to space limitations, this 
article will focus on cases where gravitational lensing is caused primarily by a 
galaxy-sized deflector (or lens). 

The first strong gravitational lens was discovered more th an thirty years ago, 



Blandford Naravan 19921 . 



decades after the phenomenon was predicted theoretically (see[ 
and references therein). However, in the past decade there has been a dramatic 
increase in the i iumber of known lenses and in the quality of the data. At the time 
of the review bv lBlandford fc: Naravan (19921 ). the 11 "secure" known galaxy-scale 



lenses could all be listed in a page and discussed individually. At the time of this 
writing, the number of known galaxy-scale lens systems is approximately 200, 
most of which have been discovered as part of large dedicated surveys with well 
defined selection functions. This breakthrough has completed the transformation 
of gravitational lensing from an interesting and elegant curiosity to a powerful 
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tool of general interest and statistical power. 

Three properties make strong gravitational lensing a most useful tool to mea- 
sure and understand the universe. Firstly, strong lensing observables - such as 
relative positions, flux ratios, and time delays between multiple images - depend 
on the gravitational potential of the foreground galaxy (lens or deflector) and its 
derivatives. Secondly, the lensing observables also depend on the overall geom- 
etry of the universe via angular diameter distances between observer, deflector, 
and source. Thirdly, the background source often appears magnified to the ob- 
server, sometimes by more than an order of magnitude. As a result, gravitational 
lensing can be used to address three major astrophysical issues: i) understanding 
the spatial distribution of mass at kpc and sub-kpc scale where baryons and DM 
interact to shape galaxies as we see them; ii) determining the overall geometry, 
content, and kinematics of the universe; iii) studying galaxies, black holes, and 
active nuclei that are too small or too faint to be resolved or detected with current 
instrumentation . 

The topic of strong lensing by galaxies is too vast to be reviewed entirely in 
a single Annual Review Article. This article is meant to provide an overview of 
a selection of the most compelling and promising astrophysical applications of 
strong gravitational lensing at the time of this writing. The main focus is on 
recent results (after ~2005). For each application, I discuss the context, recent 
achievements, and future prospects. Of course, lensing is only one of the tools 
of the astronomers' trade. When needed, I discuss scientific results that rely 
on strong lensing in combination with other techniques. For every astrophysical 
problem, I also present a critical discussion of whether strong gravitational lensing 
is competitive with alternative tools. 
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Excellent reviews and monographs are available to the interested reader for 

more details, different points of view, history of strong lensing, and a complete 
list of pre-2005 references. The Saas Fee Lectures by Schneider, Kochaneck and 
Wambsganss (2006) provide a comprehensive and pedagogical treatment of lens- 
ing fundamentals, theory and observa tions until 2 006. Additional information can 



The classic monograph by 



be found in the review bv lFalco f2005ll and that bvlCourbin Saha k Schechter f2002l') 



Schneider. Ehlers fc Falco (19921 ) and that by 



Fetters. Levine &: Wambsganss (200ll ) 



are essential references for strong gravitational lensing theory. 

This review is organized as follows. First, for convenience of the reader and 
to fix the notation and terminology, in § [2] I give a very brief summary of strong 
lensing theory. Then, in § [3l I present an overview of the current observational 
landscape. The following four sections cover the main astrophysical applications 
of gravitational lensing: "The mass structure of galaxies" (§I1D, "Substructure in 
galaxies" (§[5]), "Cosmography" (§E]), and "Lenses as cosmic telescopes" (§[7]). 
After the four main sections, the readers left with an appetite for more results 
from strong gravitational lensing will be happy to learn about the many promising 
ongoing and future searches for more gravitational lenses described in § [8j Some 
considerations on the future of strong gravitational lensing - when the number of 
known systems should be well into the thousands - are given in § [SI 

2 Brief Theoretical Introduction 
2.1 A gravitational optics primer 

Under standard conditions of a thin lens (i.e. the size of the deflector is much 
smaller then the distances between the deflector and the observer and the de- 
flector and the source), responsible for a weak gravitational field (i.e. deflection 
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angles much smaller than unity), in an otherwise homogeneous universe, strong 

lensing by galaxies can be described as a transformation from the two-dimensional 

observed coordinates associated with a particular light ray {6 in the image plane) 

to the two-dimensional coordinates that the light ray would be observed at in the 

absence of the deflector (/3 in the source plane). 

A simple and intuitive understanding of the basic principles of strong lensing by 

galaxies can be gained by co nsidering a generalized version of Fermat's principle 



(jBlandford &: Naravan 1992I . and references therein). For a given source position 



/3, the excess time-delay surface as a function of position in the image plane is 
given by: 

where Dg, D^s are respectively the angular diameter distances between the 
observer and the deflector, the observer and the source, and the deflector and 
the source, and tp is the two-dimensional lensing potential, satisfying the two- 
dimensional Poisson Equation: 

= 2k, (2) 

where k is the surface (projected) mass density of the deflector in units of the 
critical density Sc = c^Ds/i'i'^GDdDds) (for convenience of the reader I adopt 



the same notation as 



Schneider. Kochanek &: Wambsganss 20061 ). 



According to Fermat's principle, images will form at the extrema of the time- 
delay surface, i.e. at the solutions of the so-called lens equation: 

p = e - vij = e - a (3) 

which is the desired transformation from the image plane to the source plane. 
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The scaled deflection angle a is related to the deflection angle experienced by a 

light ray a hy a = ■ The lensing geometry is illustrated in Figure [TJ Note 

that the transformation is achromatic and preserves surface brightness. 

Strong lensing occurs when Equation [3] has multiple solutions corresponding to 
multiple images. Examples of the most common configurations of strong gravi- 
tational lensing by galaxies are shown in Figure [2] and explained with an optical 
analogy in Figure [3l For a given deflector the solid angle in the source plane that 
produces multiple images is called the strong-lensing cross section. For a given 
population of deflectors, the optical depth is the fraction of the sky where distant 
sources appear to be multiply-imaged. 

The Jacobian of the transformation from the image to the source plane, gives 
the inverse magniflcation tensor, which can be written as 



'dO~ dOide 



1 - K - 71 -72 
y -72 1 - K + 71 J 



(4) 



and describes the local isotropic magnification of a source (determined by the 
convergence k, defined above) and its distortion (shear components 71, and 72). 

In the limit of a point source, the local magnification /u is given by the deter- 
minant of the magnification tensor: 

(1 - k)2 - 7^ - 72^ 

For extended sources, the observed magnification depends on the surface bright- 
ness distribution of the source as well as on the magnification matrix. 

When the determinant of the inverse magnification matrix vanishes, the mag- 
nification becomes formally infinite. The loci of formally infinite magnification in 
the image plane are called critical lines. The corresponding loci in the source plane 
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are called caustics. Compact sources located close to a causti c can be magnifie d 



by very large factors up to almost two orders of magnitude (jStark et al. 20081 ). 
although the total observed flux is always finite for astrophysical sources of finite 
angular size. 

It is convenient to define the Einstein radius. For a circular deflector it is 
the radius of the region inside which the average surface-mass density equals 
the critical density. A point source perfectly aligned with the center of a circu- 
lar mass distribution will be lensed into a circle of radius equal to the Einstein 
radius, the so-called Einstein ring (see Figure [2]). The size of the Einstein ra- 
dius depends on the enclosed mass as well as on the redshifts of deflector and 



source. 



'he definition of Einstein radius needs to be modified for non-circular de- 



flectors ([Kormann. Schneider &: Bartelmann 19941 ). Once appropriately defined 



the Einstein radius is a most useful quantity to express the lensing strength of 
an object, and it is usually very robustly d etermined via strong lens models (e.g. 



Schneider. Kochanek &: Wambsganss 20061 ). As a consequence, the mass enclosed 



in the cylinder of radius equal to the Einstein radius can be measured to within 
1-2%, including all random and systematic uncertainties. 



A final essential concept is that of mass-sheet degeneracy (jFalco. Gorenstein Sz Shapiro 19^ 
Given dimensionless observables in the image plane, such as relative position, 
shape, and flux ratios of multiple images, the solution of the lens equation is 
not unique. For every mass distribution k{0) and every surface brightness dis- 
tribution in the source plane I(/3), there is a family of solutions given by the 
transformations: 

= (1 - A) + Aai; (3^ = (3/\. (6) 
The transformation changes the predicted time-delay between multiple images 
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and the magnification as follows: 

Atx = At/A; = m/A^, (7) 

resulting in a degeneracy in inferred quantities such as intrinsic luminosity and 
size of the background source. Additional information is needed to break this 
degeneracy, such as the intrinsi c luminosity or size of the lensed source (as in 



the case of lensed supernovae la lKolatt &: Bartelmann 19981 ). the actual mass of 



the deflector (as measured for example with stellar kinematics), or the measured 
time-delays between multiple images within the context of fixed cosmology. Al- 
ternatively, the mass-sheet degeneracy can be broken in the context of a model, 
for example by assuming that the surface mass density of the deflector goes to 
zero at large radii (thus A = 1). In practice, this is not always possible because 
mass structure along the line of sight - associated or not with the main deflector - 
can act effectively as a "sheet" of mass with external convergence Hcxt ■ Breaking 
the mass-sheet degeneracy is essential for a number of strong lensing applications, 
as we will see in Section [6l 

2.2 Modeling galaxies: macro, milli, and micro lensing 

It is useful to define three regimes to describe the lensing properties of the compo- 
nents of galaxies, corresponding to the typical scale of associated Einstein radii, 
as summarized in Figure HI 

2.2.1 Macrolensing On the coarsest scale, corresponding to Einstein radii 
of the order of arcseconds, the overall mass distribution of the lensing galaxy is 
responsible for the main features of the multiple images, such as image sepa- 
ration and multiplicity. In terms of physical components of an isolated galaxy, 
macrolensing can be thought of as the combined lensing properties of the DM 
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halo, the bulge, and the disk. A simple model that reproduces image positions, 

multiplicity, and fluxes is sometimes referred to as the macro model, and is gen- 
erally sufficient to infer quantities such as projected mass inside the Einstein ra- 
dius and overall ellipticity and orientation of the mass distribution. The simplest 
model that is found to provide a qualitatively g ood description of the macroscopic 



features of strong lensing by galaxies is the SIE (jKormann. Schneider Sz Bartelmann 19941 ) 
an elliptical generalization of the SIS. The three-dimensional mass density profile 
of the SIS is given by 



^sis 



^ 27rGr2 ' 
and the Einstein radius is given by 



Note that for early-type lens galaxi es asTS is found to b e approximately equal to 



the stellar velocity dispersion (e.g.. 



Bolton et al. 20n8al 'l 



An example lens model is shown in Figure [5j This system consists of a fore- 
ground elliptical galaxy lensing a background galaxy, well described by an ellipti- 
cal Gaussian surface brightness distribution in the source plane. A SIE mass 
model is found to be sufficient to reproduce accurately the observed surface 
brightness distribution in the image plane. For an SIE mass model, two curves 
(outlined in white in the figure) separate regions of different multiplicity in the 
source plane. Sources outside the outer curve (known as cut) are singly imaged, 
sources in between the cut and the inner caustic curve produce two visible images 
(plus a third infinitely demagnified central image), and sources inside the inner 
caustic produce four visible images (plus a fifth infinitely demagnified central 
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image). In this case, the detectable part of the extended source crosses the inner 

caustic so that it appears partly doubly-imaged and partly quadruply imaged. 
Due to the good alignment of the deflector and source, the image forms an almost 
perfect Einstein ring. An alternative light traces mass (LTM) model (i.e. the sur- 
face mass density of the deflector is obtained by multiplying its surface brightness 
by a mass-to- light ratio, allowed to be a free-parameter) is also shown. In this 
case, the LTM model is almost indistinguishable from an SIE model, because 
strong lensing is sensitive only to the mass enclosed by the Einstein radius, to 
first order. In general, LTM models can be excluded when considering extended 
sources because they fail to reproduce the detailed geometry, the radial behavior 
in particular. LTM models can also be excluded on the basis of a number of other 
considerations, as reviewed in § HI 

We just discussed an example of a simply-parametrized gravitational lens macro 
model, where both the source surface brightness and the mass distribution of the 
deflector are described by astrophysically motivated models with a small number 
of parameters. This kind of model is generally capable of reproducing all the 
macroscopic features while delivering robust estimates of the most important 
quantities for the deflector (e.g. total mass ellipticity and orientation) and the 
source (e.g. intrinsic size and luminosity). For these reasons, simply parametrized 
models are often all one needs in interpreting lensing data. 

However, some applications require more sophisticated lens models, capable of 
extracting more detailed information. In recent years, the increase in number of 
known lenses has been paralleled by ever more sophisticated lens modeling tools. 
A full description of advanced lens models is beyond the scope of this review. 
However, I list a few examples to point the interested reader towards the tech- 
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nical literature. A number of groups have developed "grid-based" models (also 

known - incorrectly - as non-parametric models; pixel values are parameters like 

any other), where the potential (or surface mass density) of the deflector and/or 

the surface brightness of the source are described by a set of pixels on regu- 



due to noise (e.g. 


Brewer &: Lewis 20061. 


iDve & Warren 2005. 


KooDmans 2005 




Suvu et al. 2006 


Treu & Koopmans 2004 


Veeetti & Koopmans 2009a 


Warren k Dve 2003 



An alternative hybrid approach consists of using large numbers of simply-parametrized 
models to strike a balance between flexibility and jj rior information on the shape 
and surface brightness of galaxies jnarshall 2005! ) . Bayesian statistics has be- 
come the standard statistical framework for advanced models, allowing for rig- 
orous analysis of the uncertainties in highly dimensional spaces as well as quan- 
titative model sele ction. Heuristic pixellat ed approaches have also been adopted 



with some success (ISaha fc Williams 20041 ). and recently been cast in a Ba yesian 
framework to improve the understanding of the uncertainties (jColes 20081 ). 



2.2.2 MiLLiLENSiNG On an intermediate angular scale are the lensing effects 



introduced by substructure, both luminous and dark. Typically, a lens galaxy wi 



have some satellites, like the dwarf satellites of the Milky Way (jKravtsov 20 id . 
and references therein). The mass associated with the satellites introduces per- 
turbations in an otherwise smooth potential. These perturbations can be detected 
relative to a smooth model using accurate measurements of flux ratios, relative 
position, and time delays between multiple images. This regime is sometimes re- 
ferred to as millilensing - due to the characteristic milli-arcsecond Einstein radii 
expected for dwarf satellites of massive galaxies. However, the phenomenon could 



span several orders of magnitude, depending on the mass function of satellites 
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and their spatial distribution (e.g. 



Kravtsov20ia). 



2.2.3 MiCROLENSiNG Finally, on the smallest angular scale, galaxies are 
made of stars. The Einstein radius of a solar mass star at a cosmological distance 
is of the order of micro arcseconds, hence the name cosmological microlensing. 
The average projected separation of stars in distant galaxies is small compared 
to the typical Einstein radii, and thus every background source effectively expe- 
riences cosmological microlensing. As in the case of galactic microlensing, the 
resolution of current instruments is insufficient to detect cosmological microlens- 
ing via astrometric effects. However, if the angular size of the background source 
is smaller or comparable to the typical stellar Einstein radius, cosmological mi- 
crolensing can be detected by its effect on the observed flux. In contrast, if the 
source is much larger than the typical stellar Einstein radius, the total magnifica- 
tion will be effectively averaged over a large portion of the magnification pattern 
and therefore be similar to that expected for a smooth model. The relative mo- 
tion of stars with respect to the background source and center of mass of the 
deflector are sufficiently fast to modify the magnification pattern over timescales 
of just a few years, as illustrated in Figure [6j 

As we will see in the rest of this article, all three regimes can be used to infer 
unique information on the distribution of mass in (deflector) galaxies, and on 
the surface brightness distribution of distant (lensed) galaxies and active galactic 
nuclei with sensitivity and resolution beyond those attainable without the aid of 
gravitational lensing. 
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3.1 Present-day samples and challenges 



Approximately 200 examples of strong gravitational lensing by galaxies have been 
discovered to date. A number of different strategies have been followed. The two 
most common strategies start from a list of potential sources or potential deflec- 
tors and use additional information to identify the (small) subset of strong gravi- 
tational lensing events. Other promising approache s include searching for gravita- 



tional lensing morphologies in high resolution data ( Marshall et al. 20091. and ref- 
erences therein) and exploiting variability in time domain data (jKochanek et al. 2006al ) 
The current state of the art is illustrated in Figure [71 which shows the redshift 
distribution of the lenses discovered by the four largest surveys to date. The first 
two are source-based surveys, the third is a deflector-based survey, and the fourth 
one is a lensing morphology survey. 

The Cosmic Lens All-Sky Survey (CLASS) is based on radio imaging. They 
discovered 22 multiply-imaged active nuclei, including a subset of 13 systeni s 



which are known as the statistically well-defined sample (jBrowne et al. 20031 ). 
Source and deflector redshifts are available for 11 and 17 systems, respectively 
(C.D.Fassnacht, 2009 priv. comm). The SDSS Quasar Lens Search (SQLS), iden- 



tified 28 galaxy-scale multiply-imaged quasars using SDSS multicolor 



data to sift through the spectroscopic quasar sample (jOguri et al. 20061 . 



iinagmg 



2008|) 



All source redshifts are availa ble, while deflector redshifts are available for 15 



systems. The SLACS Survey (jBolton et al. 20061 ) is an optical survey based on 



spectroscopic preselection from SDSS data and imaging confirmation with HST. 
SLACS discovered 85 galaxies acting as strong lenses (plus an additional 13 prob- 
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able lenses; 
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Auger et al. 2009i ). Source and deflector redshifts are available for all 



systems. Finall y, twenty secure galaxy-scale len s systems were discovered by vi- 



sual inspection (jFaure et al. 20081 . 



Jackson 20081 ) of the HST images taken as part 



of the COSMOS Surv ey. Source and deflecto r redshifts are available for 3 and 13 



systems, respectively (jLagattuta et al. 20091 ). 



The compilation is not complete, due to the difficulty of keeping track of the 
ever gr owing number of lenses discovered serendi^ itously or by ongoing concerted 
efforts fjCabanac et al. 20071 . IMarshall et al. 20091 ) that still lack confirmation and 
spectroscopic redshifts (a useful resource to find data for lenses from a variety of 
sources is the online database of strong gravitational lenses CASTLES at URL 



http://www.cfa.harvard.edu/castles ). However, the compilation gives a good idea 



of the observational landscape and of the two main limitations of current samples. 
Firstly, most new lenses have been found at z < 0.4, which is a very favorable 
regime for detailed follow-up, but limits the look-back time baseline for evolution- 
ary studies and the spatial scales probed by lensing. Secondly, many gravitational 
lens systems still lack source or deflector redshifts. 

It is customary to classify strong lenses as galaxy-galaxy lenses (e.g. Figs [2] and 
[5]), and galaxy-QSO lenses (e.g. Figure [6]), depending on whether an active galac- 
tic nucleus is present in the backgroun d source. Galaxy-QSO lenses are in ore rare 



on the sky than galaxy-galaxy lenses (jMarshall. Blandford &: Sako 20051 ). How- 



ever, they can be found efficiently by exploiting their radio emission and the 
variability of the point source. Furthermore, the compact point source enables 
studies of the granularity of the lens galaxy (from microlensing) , and of cos- 
mography and lens galaxy structure (from direct measurements of time delays 
between images). Galaxy-galaxy lenses are typically more suited for the study 
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of the lens galaxy itself, because its emission is not overwhelmed by the multiple 

images of the background source. Furthermore, the extended surface brightness 

of the source provides detailed information on the gravitational potential of the 

deflector. 

It is observationally challenging to extract the wealth of information available 
from strong lensing systems. First and foremost, subarcsecond angular resolution 
is key to identifying and characterizing strong lensing systems. Radio or opti- 
cal/near infrared observations from space (and recently from the ground with 
adaptive optics) have been essential for the progress of the field. Secondly, both 
source and deflector redshifts are needed to transform angular quantities into 
masses and lengths. Especially for the so urce redshift, lon g exposure on the 



Qfek et al. 20061 ). Success is not as- 



largest telescopes are typically required (e.g. 
sured, and in many cases one must rely on photometric redshifts, which are also 
challenging because light from the foreground deflector complicates photometry 
of the background source. Third, micr olensing and variab ility depend critically 



'his makes X-ray (e.g.. 



Poolev et al. 20091 ) and mid-infrared ob- 



on source size. 

servations (e.g., lAgol et al. 20091 ) ~ probing sources that are much smaller and 
much larger than the scale of microlensing, respectively ~ particularly useful, even 
with limited spatial resolution. Fourth, time delays and microlensing studies re- 
quire intensive monitoring campaigns, with all the associated logistical challenges. 
Lastly, depending on the application, ancillary data such as velocity dispersion 
or information on the local large scale structures are typically needed to break 
degeneracies and control systematic errors. 
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3.2 Selection 

Strong lensing is a very rare phenomenon. Wi th present technology only ^1/100 



galaxies can be detected as strong lenses (jMarshall. Blandford fc Sako 20051 ). 



Similarly, the optical depth is of order 10 -10 , i.e. < 1/1000 high-redshif t 



Browne et al. 20031 ) 



sources in the sky have detectable multiple images (e.g. 
Both numbers depend strongly on the depth of the observations. Thus, in order 
to generalize the results obtained from this technique to the overall population 
of deflectors and sources, and for applications of strong lensing to cosmography, 
it is essential to understand the selection function very well. 

To first order, strong lensing galaxies can be described as selected by veloc- 
ity dispersion. Most galaxy-scale strong gravitational lenses discovered to date 
are massive elliptical galaxies with velocity dispersions in the range 200-300 
kms~^. This well-understood selection function arises from the rapid increase 
in the strong lensing cross section with mass (oc o"^ for an SIS), and from the 
rapi d decline of the velocity dispersion functi on of galaxies above 300 kms~^ 



(see 



Schneider. Kochanek fc Wambsganss 2006t . for a comprehensive discussion). 



As an example, the average stellar velocity dispersion of the SLACS sample is 
248 kms~^ with a r.m.s. scatter of 46 kms ^. The velocity dispersion selection 
is also responsible for the adverse selection against late-type galaxies. Approxi- 
mately 80% of the SLACS deflectors are pu re ellipticals, 10% are lenticulars and 



10% are spirals, mostly bulge dominated (j Auger et al. 20091 ). Identifying and 
studying galaxies with a < 200 kms^^ acting as strong gravitational lenses is 
possible with sufficiently large surveys and represents an exciting frontier for the 
next decade. However, this is an observationally challenging problem because 
the image separation drops quickly below 0'.'3 — 0'.'4, the current practical limit 
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for detection with HST and the Very Large Array (VLA). Furthermore, once the 

resolution drops significantly below the typical arcsecond size of distant galaxies, 
disentangling light from the deflector and background source becomes increas- 
ingly difficult, particularly at optical/infrared wavelengths. 

The lensed sources are to first order fiux and surface brightness selected. This 
translates into a complex selection function in terms of the intrinsic properties of 
the source population because of the magnification effects of lensing. It is easier 
to understand the effect for point source surveys, such as CLASS and SQLS. Due 
to lensing magnification, sources that are fainter than the survey flux limit will 
enter the sample. However, magnification also reduces the solid angle actually 
surveyed. Therefore, the number of strong lensing events depends critically on 
the dependency of the surface density of sources on the observed flux. This effect 
is known as magnification bias. For extended sources, observed magnification will 
also depend on surface brightness and size of the source, generally being larger 
for more compact sources. The redshift distribution of the lensed sources will in 
general be different than that for a non-lensed population selected to the same 
apparent magnitude limit. 

Other more subtle selection effects are also at work. Factors that may affect 
the strong lensing cross-section of a galaxy include elongation along the line of 
sight, flattening of the projected mass distribution, concentration of the mass 
distribution (e.g., the slope of the mass density profile at fixed virial mass), 
overdensity of the local environment, and abundance of small scale structure in 
the plane of the deflector or along the line of sight. Factors that may affect the 
probability of a source being identified as multiply-imaged include extinction from 
the foreground lens galaxy, configuration of the multiple images (in particular 
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image separation and flux ratios), time variability, and presence of emission lines 

and hence properties of the stellar populations or existence of an active nucleus. 

Three complementary strategies have been followed to quantify selection ef- 
fects. One strategy consists of starting from a realistic cosmological model and 



Mandelbaum. van de Ven Keeton 20091 . 



simulating the selection process from first principles (e.g. 
and references therein) . This is the most direct way to compare observations with 
theoretical models of galaxy formation. The challenge of this approach is that 
lensing selection depends on the details of the mass and surface brightness distri- 
butions on scales much smaller than a galaxy. Unfortunately, realistic simulations 
of the universe on this scale are beyond our current capabilities. Therefore, one 
needs to rely on DM-only simulations and approximate the effects of baryons, 
with all associated uncertainties. A second strategy consists of comparing sam- 
ples of lens galaxies with control samples of non-lens galaxies. This approach was 
used with the SLAGS sample to show that - once velocity dispersion and red- 
shifts are matched - lens galaxies are indistinguishable within the uncertainties 
from twin galaxies selected from SDSS in terms of their size, surface brightness. 



lumin osity, location on t 



ment (jAuger et al. 20091 . 



re fundamental pla,ne, stellar mass, and lo cal environ- 



Bolton et al. 2008a 



Treu et al. 2009 



2009 1. This find- 



ing implies that the results from the SLAGS Survey can be applied to the overall 
population of velocity dispersion selected early-type galaxies. The strength of 
this method is its ability to take into account real selection functions with all the 
inherent complexity. This guarantees that one compares apples with apples, but 
does not solve the problem of comparing with theoretical models. A "hybrid" 
approach consists of constructing simple models starting from empirically-based 
information on the deflector and source populations, and combining it with lens- 
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ing theory to compute the relevant selection function. This approach is extremely 

useful for developing an intu ition for the process and compute approximate cor- 
rection factors. For example, 



1 

Oguri (2007i ) was able to explain the observed ratio 



of quadruply imaged to doubly imaged quasars in the CLASS sample in terms 
of magnification bias. The challenge for this approach is including a sufficiently 
accurate description of the physics and details of the observations to infer quan- 
titatively correct answers. 

4 The Mass Structure of Galaxies 

The standard cosmological model, based on CDM and dark energy reproduces 
very well the observ ed structure of the universe on supergalactic scales (e.g.. 



Komatsu et al. 2009i . and references therein). At galaxy scales, DM and baryons 



interact to produce the observed variety of galaxy properties. The situation is 
not so clear at small sub-galactic scales, wher e potential conflic ts between theory 



Ellis k Silk 20091 ). Understanding 



and observations have been suggested (e.g. 
the interplay between DM and baryons is crucial to make progress in developing 
and testing theories of galaxy formation at these scales. Gravitational lensing, 
by providing direct and precise measurements of mass at galactic and subgalactic 
scales, is a fundamental tool for answering a number of questions with profound 
implications on the existence and nature of DM. Do galaxies reside in DM halos? 
How do the properties of galaxies depend on those of their DM halos? Are DM 
density profiles universal as predicted by simulations? These are the topics of 
this section. 



Strong Lensing by Galaxies 21 
4.1 Luminous and dark matter in early-type galaxies 

4.1.1 Do EARLY-TYPE GALAXIES LIVE IN DARK MATTER HALOS? It is com- 
monly believed that all galaxies live in DM halos. However, in the case of early- 
type galaxies, observational evidence is hard to obtain. The difficulty arises 
mostly from the paucity of mass tracers at radii much larger than the effective 
radius Re - where DM dominates - and from the degeneracies inherent in in- 
terpreting projected data in terms of a three-dimensional mass distribution for 
pressure supported systems. Chief among these degeneracies is that between 
the total mass density profile and the anisotropy of th e pressure tensor ("mass- 



anisotropy" degeneracy, e.g., 



Treu k Koopmans 2002al l. 



Much progress in detecting DM halos has been achieved by studying the 
kine matics of stars, globular clusters, and cold and hot gas in nearby systems 



(e.g., 



Bertin fc Stiavelh 1993 



Humphrey et al. 20061 ). This type of study shows 



that DM halos are generally required to explain the dynamics of massive early- 
type galaxies. Weak-lensing has been used to demonstrate the existence and 
to characterize the outer regions of DM halos for statistic al samples of early- 



type galaxies out to intermediate redshift s {z ~ 0.5, e.g. 



Hoekstra et al. 2005 



Lagattuta et al. 20091 ) 



Gavazzi et al. 20071 . 



Strong lensing observations demonstrate the existence of DM halos around 
individual massive early-type galaxies out to z ~ 1 beyond any reasonable doubt, 
both by themselves and in combination with other techniques (for early-type 
galaxies with a < 200 kms~^ the case is much less conclusive; future sample of 
low-mass deflectors may be needed to clarify matters). One argument is that the 
amount of mass inside the Einstein radius exceeds the stellar mass M,, . This latter 
quantity can be constrained in many ways. Assuming an IMF, stellar population 
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synthesis (SPS) models applied to photometric or spectroscopic data yield M* 

with an uncert ainty of 0.1-0.2 dex. Al t ernatively, local dyna mical studies of early- 



type galaxies (jCappellari et al. 20061 . iGerhard et al. 20011 ) constrain the stellar 



mass-to-light ratio at present time, which can then be evolved back in time either 
using the measured evolution of the fundamental plane or other measurem ent of 



the star formation history (e.g.. iKochanek 19951 . iTreu &: Koopmans 2004 ) 



A particularly powerful combination for detecting DM halos is to use stel- 
lar kinematics of the lens galaxy to provide information on the distribution 
of mass in the high surface brightness regions well within the effective radius, 



and to use strong 



e nsing to help remove th e mass-anisotropy degeneracy (e.g. 



Barnabe et al. 2009 



Treu &: Koopmans 2004 ). A third method relies on assum- 



ing scaling relations to analyze lenses across a sample and reconstruct the mass 
density profile for the ensemble, which turns out to be more extended than ex- 



pected if mass followed li^ ht and therefore consistent with DM (iBolton et al. 2008b 
Rusin Kochanek 2005l ). A fourth method exploits microlensing statistics to 
demonstrate that point masses (i.e. stars) cannot contribute the totality of the 



Poolev et al. 20091 ^. 



surface mass density at the location of the multiple images (e.g. 
A fifth method consists of measuring time delays between multiple images, de- 
termining angular-diameter distances from independent cosmographic probes to 
i nfer the behavior of the mass density profile at the location of the multiple images 



(IKochanek et al. 2006bl ^. 



4.1.2 What is the relative spatial distribution of luminous and 
DARK MATTER? The efficiency with which baryons condense inside halos to 
form stars, and their effect on the underlying DM distribution, depend on the 
interplay between cooling and heating (e.g. from star formation and nuclear 
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activity). Lensing can help us understand these process by providing precise 

measurements of the fraction of total mass in the form of DM (foM) within a 

fixed proje cted radius, typically exp ressed in terms of fraction of the effective 



radius (e.g. 



Jiang fc Kochanek 200?! ) 



Qbservatio nally, fpM is found to be non- negligible already at the effective ra- 



dius (25±6 



at five effective radii 



KooDinans et al. 20061) and incr easing towards larger radii (70ibl0% 



Treu &: Koopmans 2004 ) . Consis tent results are obtaine d by 



a number of independent non-lensing techniques (e.g. iCappellari et al. 20061 ). In 
addition, foM within a fixed fraction of the effective radius is found to increase 
with galaxy stellar mass and velocity dispersion. For example, by comparing 
lensing masses with those inferred from SPS modeling of multicolor data, fDM 
inside the cylinder of projected radius equal to the Einstein radius increases from 
~ 25% to ~ 75% in the range of velocity dispersion cr=200-35 kms~^, or equiv- 



alently in the range of stellar mass between 10^^ and 10^^ (lAueer et al. 20091 . 



I 1 

Figure [8]). These numbers are based on a iSalpeter f 19551) IMF and are consis - 



tent with t hose inferred by local dynamical studies (e.g 



Cappellari et al. 20061 ). 



Adopting a iChabrier (20031) IM F changes the overall normalization, but not the 



global trend ( 



Auger et al. 20091 . see, however, Grillo et al. 2009 for a contrasting 



view . 



Strong lensing stud ies also explain the origin of th e so-called tilt of the fun 



damental plane (e.g., 



Ciotti. Lanzoni &: Renzini 19961 . hereafter FP), the tight 



correlation between effective radius, effective surface brightness and stellar ve- 
locity dispersion observed for early-type galaxies. By introducing a dimensional 
mass variable Mdim = o"^Rc/G, the FP can be cast in terms of an increasing effec- 
tive mass-to- light ratio with effective mass (the tilt'). Exploiting strong lensing, a 
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somewhat tighter mass plane (MP iBolton et al. 2008bl ) relation can be obtained 
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by replacing surface brightness with total surface mass. The MP is not tilted, 
implying that the tilt of the FP stems from an increase in foM with mass, and 
not in a systematic change, e.g., of the virial coefficient that connects Mdim to 
total mass. 

4.1.3 Mass density profiles and the bulge-halo conspiracy An- 
other quantity of interest is the average logarithmic slope of the three-dimensional 
total mass density profile d log ptot log r = —7'. An isothermal mass model has 
7' = 2. The total mass density profile for a spherical model is often expressed in 
terms of the equivalent circular velocity 



^^B<ll, (10) 

r 

which facilitates comparison with the literature on spiral galaxies and on numer- 
ical simulations. For a spherical power-law density profile, 7' is simply related to 
the slope of the rotation curve by the relation dlogVc/dlogr = (2 — j')/2. For 
this reason, an isothermal profile is sometimes referred to as a flat rotation curve. 

The basic result on this topic is that 7' w 2, i.e. early-type lens galaxies 
have approximately isothermal mass density profiles, or close-to-flat equivalent 
rotation curves. This has been known since at l east the early nineties, both 



on the basis of lensing studies (e.g. 



1 r 



Kochanek 19951 ) and on local kinematics (e.g. 



Bertin &: Stiavelli 19931 . iGerhard et al. 200ll . and references therein). However, in 
order to understand the mass structure of galaxies with sufficient level of precision 
to constrain formation models, we need to ask more detailed questions. What 
is the average 7' and its intrinsic scatter for the overall population of early-type 
galaxies? How does 7' depend on the galactic radius or other global properties? 
Does it depend on the environment, as expected if halos were tidally truncated? 
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Does 7' evolve with redshift? In addition, as we will see in § [6l determining the 

mass profiles of lens galaxies to high accuracy is essential for many applications 

to cosmography. 

In the past few years, the large number of lenses discovered and the high 
level of precision attainable with lensing has enabled substantial breakthroughs. 
Joint lensing and dynamical studies of the SLACS sam ple have shown that 7^ 



2.08 lb 0.02 with an intrinsic scatter of less than 10% (jKoopmans et al. 2009bl ) 



This result is valid in the sense of an average slope inside one effective radius 
or less, the typical size of the Einstein radius of SLACS lenses. For higher red- 
shift deflectors, Einstein radii are typically larger than the effective radius and 
reach out to 5 Rg. Although the high redshift samples with measured velocity 
dispersions are small, they seem to sug gest a somewhat larger intrinsic scatter 



around 7' = 2 (jTreu Koopmans 2004 ). No significant dependency on galactic 



radius, global galaxy param eter, or redshift has been found so far based on lens- 
ing and dynamical analysis (jKoopmans et al. 2009bi ). The small scatter around 



7' = 2 is remarkable, considering that neither the DM halo, nor the stellar mass 
are well described by a simple power-law profile. Nevertheless, the two com- 
ponents add up to an isothermal profile (Fig. [9]). This effect is similar to the 
disk-halo conspiracy respoii sible for the flat rotation curves of spiral galaxies 



(jvan Albada &: Sancisi 19861 ) . and it is therefore been dubbed the 'bulge-halo 
conspiracy'. Detailed dynamical studies of the twodimensional velocity field of 
deflector galaxies in conjur iction with strong g ravitational lensing confirm this 



picture to higher accuracy ([Barnabe et al. 20091 ). 



Similar and consistent results can be obtained directly from gravi tational lens 
models, both for lensed sources covering a significant radial range (e.g. 



Dve fc Warren 20051 1 
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or when a gravitational time del ay has been measured a nd the cosmology is fixed 



by independent measurements ([Kochanek et al. 2006bl ). An interesting case is 



that of the system SDSSJ0946+1006 where the presence of two multiply- imaged 
sources at different redshifts constrains the projected mass density slope to be 
7' = 2.00 lb 0.03, based purely on lens modeling (Figure [T0|) . The lack of cen- 
tral images also constrains the slope of the total density profile to be steep (e.g., 
7' = 2) in the central regions of deflectors. It should be noted that lensing is 
mostly sensitive to the projected mass density slope at the location of the im- 
ages, rather than the average inside the images. Therefore, a direct comparison 
with the lensing and dynamical results is only valid to the extent that a pure 
power-law profile is a good model for the data. 

4.1.4 Are dark matter density profiles universal? Cosmological 
numerical sim ulations predict that DM densit y profiles should be almost universal 



in their form (jNavarro. Frenk fc White 19971 . hereafter NFW) . Simulated profiles 



are characterized by an inner slope dlog pDAi/dlogr = —7 —1. At the scales 
of spiral galaxies, low surface brightness galaxies, and clusters of galaxies, it has 
been shown that in a nur nber of systems the observed profiles are shallower than 



predicted (i.e. 7 < 1, e.g., 



Salucci et al. 2007 



Sand et al. 20081 ) . The discrepancy 



suggests that either the DM component or the effects of baryons on the underlying 
halos are poorly understood. 

In early-type galaxies the inner regions are completely dominated by stellar 
mass, making them particularly interesting systems for understanding the in- 
terplay between baryons and DM. Unfortunately, the dominance of baryons also 
makes the measurement more challenging. A joint lensing and dynamical analysis 
of 5 high-z lenses shows that 7 is consistent with unity, albeit with large errors. 
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and shallower slopes cannot be excluded (jTreu &: Koopmans 2004 ). Improving 

the measurement will require larger samples of objects with good quality data 

and further constraints on the stellar mass-to-light ratio. 

Alternatively, by imposing 7 = 1 one can infer an absolute normalization of the 

stellar mass component, an d thus constrain the IMF of massive early- type galax- 



ies to be close to Salpeter (jGrillo et al. 2009 . 



Treu et al. 2010l ). A joint lensing. 



dynamical, and stellar population analysis of the SLAGS sample shows that mas- 
sive early- type ga laxies cannot have both a universal DM halo and universal IMF 



(jTreu et al. 2010l ): either the inner slope of the DM halo or the normalization of 
the IMF have to increase with deflector velocity dispersion. 
4.1.5 Implications for early-type galaxy formation Massive early- 
type galaxies are simple dynamical systems with simple stellar populations. Yet, 
their formation and evolution is still far from being well understood (for a com- 
prehensive review see S.M. Faber, this volume). The standard CDM model pos- 
tulates their formation via major mergers, but this is hard to reconcile with their 
uniformly old stellar populations - unless there is some fine-tuned f eedback mech- 
anism that prevents star formation in the high mass systems (see lRenzini 20061 . 



for a recent review) - and with the slow observed evolution of their stellar mass 
function since z ~ 1. Recently, collisionless mergers not involving gas and star 
formation (and therefore "dry") have become increasingly popular as a possi- 
ble mechanism of growth (Faber, this volume). Furthermore, dry mergers can 
grow galaxies in size faster than in velocity dispersion. Therefore they have 
been suggested as a possible mechanism for the evolution of ultradense massive 
galaxies at high redshift into the more diffuse ones found in the local universe 



(jvan der Wei et al. 200% . 
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Strong lensing studies give us some direct information on the connection be- 
tween baryons and DM, and therefore offer us new insights into this problem. 
The (non-evolving) isothermality of the total mass density profile requires an 
early dissipative phase, to steepen the NFW profiles predicted in CDM-only sim- 
ulations. Alternatively, an initial collapse associated with incomplete violent 
relaxation could have established the isothermality of the inner profiles. Either 
phenomenon must have occurred well before z ~ 1. After the initial formation, 
further growth by dry mergers preserv e s the isothermal profile and tightness of 



the mass plane (jKoopmans et al. 20061 : 



Nipoti. Treu Bolton 2003 ). However 



dry mergers do not preserve the tight c orrelations between size and total mass 



and velocity dispersion and total mass (jNipoti. Treu &: Bolton 2003 ). The ob- 



served tightness of the correlation limits the growth by dry mergers to have been 
at most a factor of two since z ~ 2, unless there is a large degree of fine tuning 
between orbital parameters of the merger and location in the size-mass-velocity 
dispersion space. Therefore, it seems most likely that the majority of the mass 
assembly must have occurred during the initial dissipative phase associated with 
the dominant episode of star formation. 

The other main strong lensing result, i.e. the correlation between DM frac- 
tion and velocity dispersion (stellar mass), provides us with anot her piece of the 



puzzle. Dry mergers increase fdm (|Nipoti. Treu &: Bolton 2003 ). thus creating 



part of the trend. However, dry-mergers cannot explain the whole trend, which 
must be largely established early-on through other means. A scenario where the 
time since major initial collapse increases with present-day mass could explain 
the trend in terms of the evolution of the density of the universe with cosmic time 



([Thomas et al. 20091 ). The correlation between present day mass and epoch of 
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major mass assembly could also help explain the correlations be tween present da y 



mass, age, and chemical composition of the stellar populations (jTreu et al. 20051 ). 

It should be noted that the conclusions above hold only for the most massive 
early-type galaxies. At lower masses, evolution is certainly more recent and other 
secular or environmentally driven mechanisms cou ld be responsible for forming 



early-type galaxies (e.g. 



Bundv. Treu fc Ellis 2007l 'l 



4.2 Luminous and dark matter in spiral galaxies 



Massive DM halos around loc al spiral galaxies are readi ly detected from the 



gas kinematics at large radii ([van Albada fc Sancisi 19861 ). The total gravita- 



tional potential can be reconstructed accurately from the observed velocity field. 
However, decomposing the total mass distribution into its baryonic and dark 
components for individual galaxies is still an unsolved problem, largely because 
the stellar mass-to-light ratio is uncertain by a factor of ~2-3 for young and 
dusty stellar populations. In the distant universe, the problem is compounded by 
observational difficulties: HI becomes prohibitively expensive to detect; optical 
rotation curves can be measured out to z ~ 1 but are limited by cosmological 
surface brightness dimming as well as angular resolution. One approach consists 
of assuming that the b aryonic component is maxim ally important, the so-called 



maximum-disk ansatz (jvan Albada &: Sancisi 19861 ). However, it is not clear that 



disks are indeed maximal . Indeed, submaxima 
variety of arguments (e.g. 



is a dominant source of uncertainty (jBell &: de Jong 2001 



disks seem to be suggested by a 



Courteau Sz Rix 19991). even thou gh the unknown IMF 



). Understanding the 



relative mass in disks and hal os is critical to form ulate and test a robust theory 
of disk galaxy formation (e.g. iDutton et al. 20071 ). 



30 Tommaso Treu 

Gravitational lensing provides a new tool for luminous and dark matter decom- 
position in spiral galaxies. Two factors make lensing particularly useful in this 
respect. Firstly, it measure the total projected mass within a cylinder. This can 
then be combined with the enclosed mass in 3D inferred from disk kinematics 
to break the disk-halo degene racy by exploiting the different radial dependency 



of the two components (e.g., 



Mailer. Flores &: Primack 1997l ). Secondly, gravi- 



tational lensing provides azimuthal information which also helps pin down the 
relative contribution of the two, especially if they are misaligned. 

Strong lensing studies of spiral galaxies have shown encouraging results, al- 
though the impact of the conclusions is limited by the small size of current sam- 
ples. For example. 



Trott et al. (2010l ) combined lensing constraints, high reso- 
lution imaging data, and optical and radio kinematics to decompose the mass 
profile of the Einstein Cross lens galaxy into its bulge, disk, and halo components 
(see also van de venn et al. 2010, submitted to ApJ). The mass-to-light ratio 
of the bulge is very well constrained (M/L^ = 6.6 it 0.3 in solar units). Due to 
the unusually small Einstein radius of this system, the mass of the disk is less 
well constrained, although it is clearly sub-maximal, contributing 45 it 11% of 
rotational support at 2.2 scale lengths. 

The situation is changing rapidly, due to progress in strong lensing searches. 
SLACS discovered approximately 7 new bulge-dominated spiral lenses and an on- 
going search based on a similar-strategy (SWELLS; HST-GO- 11978) should find 



as many edge-on lat e type spirals. Dedicated searches (e.g.. 



Feron et al. 20091 . 



Marshall et al. 20091 ) should discover tens of new systems in the next few years. 
At variance with the smoothness of early-type galaxies, the small-scale structure 
of the surface brightness of the spiral lens due to dust and inhomogeneous stel- 



Strong Lensing by Galaxies 31 
lar populations complicates the identification and modeling of multiply-imaged 

parts of the background source. High resolution near- infrared images with adap- 
tive optics or with HST and JWST, coupled with multicolor optical data, or in 
the radio, will be essential to make progress on this front (Fig. [TT]) . 

5 Substructure in Galaxies and the "Excess Subhalos" Problem 
5.1 Background 

In the standard cosmology, DM halos host a hierarchy of sub-halos, also known 
as DM substructure. The number of subhalos above a given mass scales approx- 
imately as the total mass of the parent halo, and the logarithmic slope of the 
subhalo mass function is approximately d N/dM^„-h oc M^"^, with a = 1.9 it 0.1 



(Diemand et al. 200S 



Springel et al. 20081 ). Remarkably, the normalized distri- 



bution of substructure depends very little on the overall scale of the halo, therefore 
we would expect approximately the same abundance of satellites around clusters 
and galaxies. 

Although realistic simulations including baryons and non-gravitational effects 
have yet to be performed at this scale, it is currently believed that the sta- 
tistical properties of the substructure inferred from N-body simulations should 
be robust enou gh to allow for a direct comparison with observations (see, e.g., 



Kravtsov 201Q . and references therein). For these reasons such a comparison 



may provide one of the most stringent and direct tests of the CDM paradigm at 
subgalactic scales. 

At variance with the results of simulations, the abundance of luminous satellites 
observed around real clusters and galaxies are very different. Whereas clusters 
of galaxies host thousand of galaxies within their own DM halos, fewer satellites 
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are generally seen around galaxies. In particular, the mass function of the lumi- 
nous satellites of the Milky Way differs dramatically from that of the subhalos 
of a typical simulated halo of comparable mass. At the high mass end of the 
distribution (virial Mgub ~ 10^ Mq) the observed number of satellites is compara- 
ble, or perhaps even slightly larger, than expected. However, the mass function 
of the halos of the observed satellites is found to be much shallower than that 
predicted for subhalos, resulting in a dramatic shortfall at lower masses, below 
10^ Mq. This d iscrepancy between theory and observ ations has been known for 



over a decade (jKlvpin et al. 1999 



Moore et al. 19991 ). and has not been solved 



by the revolutionary discovery of low luminosity satellites of the Milky Way by 
SDSS, nor by advances in numerical sim ulations. An up- do-date summary of the 



current state of the problem is given by 



Kravtsov f20ld ). 



5.2 Possible solutions 



There are two classes of possible explanations for the so-called "excess sub-halos 
problem" (or "missing satellites problem" if you are a theorist). One possible 
explanation is that substructure exists, but it is dark, i.e. subhalos do not form 
enough stars to be detected. This explanation would imply that the conversion 
of baryons into stars is inefficient for small halos. It is hard to explain this 
inefficiency with the known mec hanisms of sup ernovae feedback or the effect 



of the UV ionizing background ([Kravtsov 20101 ) . Alternatively, it is possible 



that subhalos are not as abundant as predicted by numerical simulations. This 
explanation would imply a major revision of the standard CDM paradigm, either 



reducing the amplitude of fluctuat ions on the scales of sate. 



nature of DM from cold to warm ([Miranda &: Maccio 20071 ). Either explanation 



lites, or changing the 
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has far reaching imphcations. In order to be viable, the first explanation requires a 

clear improvement in our understanding of galaxy formation. In its most extreme 

version, the second explanation may require a re-thinking of the paradigm. 

Gravitational lensing provides a unique insight into this problem, since it is 
arguably the only way to detect dark substructure, measure its mass function, and 
compare it with the prediction of CDM numerical simulations. Even if advances 
in theories of galaxy formation could explain the luminosity function of Milky 
Way satellites, there would be still be a robust and falsifiable prediction of large 
numbers of darker satellites to be tested. 

If the mass function of sub-halos turns out to be different than that predicted by 
simulations a major revision of the theory would be required, possibly requiring 
warm DM, although it is not clear than that would necessarily be a compatible 
with all other constraints (see iKravtsov 20ld . and references therein) . 



5.3 Flux ratio anomalies 

The most striking and easiest to detect lensing effect of substructure is the pertur- 
bation of the magnification pattern. Since magnification depends on the second 
derivative of the potential, a small local perturbation can introduce dramatic dif- 
ferences in the observed surface brightness of the lensed source, without altering 
significantly the overall geometry of the system. For point sources, the presence 
of substructure results in ratios of the fluxes of multiple images that are signifi- 
cantly different than what would be predicted by a smooth macro model. This 
effect is often referred to as the anomalous flux ratios pheno menon, and has been 
used to infer the presence of substructure in lens galaxies 



Chiba 2002 



Dalai Kochanek 20021 . 



Mao fc Schneider 1998 



fiBradac et al. 20021. 



Metcalf fc Zhao 20021 ). 
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In an influential paper, 



1 
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Dalai fc Kochanek (2002i ) analyzed radio data for a sam- 



ple of seven quadruply imaged sources, and reported the detection of a surface 
mass fraction in the form of substructure between 0.6% and 7%. This observed 
fraction appears to be even higher than t he mass fraction in substructure at the 



Einstein radius predicted by simulations (jMao et al. 20041 . 



Xu et al. 20091 ). 



Substantial efforts have been devoted to investigate whether satellite-size ha- 
los are the most likely explanation of the observed flux ratio anomalies. Indeed, 
flux ratio anomalies could also arise from other effects such as microlensing - if 
the source is sufficiently compact - or a non-uniform interstellar medium which 
could variously affect light propagating along different paths. However, both 
contaminants are wavelength dependent, while flux ratio anomalies due to the 
substructure are achromatic. Therefore, observations at multiple wavelengths, 
especially radio, narrow emission lines, and mid-infrared, can be used to show 
that the the anomalous flux ratios are effectively due to substructures on scales 



much larger than stars (e.g . JAgol. Jones fc Blaes 2000l : 



Kochanek Dalai 2004 



Moustakas &: Metcalf 20031 ). Angular structure in th e macro model has b een sug- 



gested as a possible cause for flux ratio anomalies (jEvans fc Witt 20031 ). How- 
ever, in the cases when enough azimuthal information is available it has been 
shown that the angular s tructure of lens galaxies is fairly simple and well ap- 



proximated by an ellipse (jYoo et al. 20061 ). Elegant arguments based on the lo- 



cal curvature of the time-delay surface near the multiple images have also been 
used to show that anomalous flux ratios are indeed due to mass substructure 



(Chen 2009 



Kochanek fc Dalai 20041 ). A final source of concern is potential con- 



tamination from substruct ure along the line of sight, which could mimic the effects 



of true galactic satellites (jChen 2003 : 



Chen. Kravtsov fc Keeton 20031 ). Line-of- 
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sight contamination is most likely not the main cause of the anomalies observed 

so far. However, it is clear that line-of-sight contamination needs to be better un- 
derstood and quantified in order to extract the maximum amount of information 
from this powerful tool. 

An important questio n is whether the detected substructure is dark or lumi- 
nous. In some cases (e.g. 



Koopmans et al. 20021 : 



MacLeod. Kochanek fc Agol 20091 : 



McKean et al. 20071 ) it has been shown that mass associated with luminous satel- 



lites can explain the observed anomalies. Whether lumin ous substruc ture can ex- 



plain all the known anomalies is still a matter of debate (jChen 2009 ). On a case 
by case basis, the role of luminous satellites is difficult to quantify because they 
are hard to detect in the vicinity of the bright lensed quasars, where they would 
be most effective in introducing anomalies. In addition to high resolution HST 
or adaptive optics images, an accurate determination of the luminosity function 
and spatial distribution of luminous satellites of (non-lensing) massive galaxies 
may be a way to make progress. The challenge is to collect large enough sam- 
ples of non-lenses while carefully matching the selection process of the sample of 
lenses. It is important to stress that the detection of optical counterparts does 
not undermine the quest for substructure using gravitational lensing. Measuring 
the mass function of satellites - whether they are visible or not - is essential to 
test the CDM paradigm. Comparing the satellite mass function with their lumi- 
nosity function will only help in answering some of the questions related to the 
mechanisms which regulate star formation. 

The detection of substructure via anomalous flux ratios is an example of the 
power of gravitational lensing in measuring the distribution of mass in the uni- 
verse. However, the strong lensing studies to date suffer from two fundamental 
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limitations, which need to be overcome in order to make progress. The first Hm- 

itation is poor and uncertain statistics. Not only is the number of systems that 

can be used to study anomalous flux ratios tiny, but the selection function is 

poorly characterized. Therefore, the uncertainties are large and the results could 

be biased. The second major limitation is the limited mass sensitivity achieved 

so far, which is only sufficient to probe the upper end of the mass function of 

subhalos. 

Major improvements on both aspects are underway and significant progress is 
possible in the next few years. One key factor is the increase in the number of 
known lenses, discovered with a well defined selection algorithm, coupled with 
the increased capability for follow-up. In the next decade, we may expect tens 
of thousands of lenses to be discovered by radio and optical surveys (§ The 
other key factor is the development of advanced techniques to be applied to high 
resolution data to probe further down the mass function of subhalos, discussed 
next. 

5.4 Astrometric and time-delay anomalies 

Flux ratio anomalies is only one way to detect substructure. Subhalos affect all 
lensing observables, including deflection angles and time delays, and can therefore 
be detected as corresponding perturbations with respect to the predictions of a 



smooth model. Although these are more subtle effects, they 
to be sufficiently large to b e used to detect substructure (e.g. 



rave been shown 



Chen et al. 2007I . 



Keeton &: Moustakas 20091 ). Galaxy-galaxy lenses where the multiple images 
form an almost complete Einstein ring and are observed with high signal-to- 
noise ratio can detect individual substructures with masses as low as ~10® M.q 
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(jKoopmans 2005 



1 

Vegetti &: Koopmans 2009al ). Recent calculations bv lVegetti &: Koopmans (2009bl ) 



indicate that current samples of galaxy-galaxy lens systems such as SLACS can 
detect subhalo mass fractions as low as 0.5%, assuming the slope of the mass 
function is well known from simulations. A sample of 200 Einstein rings with 
data of comparable quality to HST should be sufficient to start constraining the 
slope of the mass function as well. The sensitivity will be further enhanced with 
advances in resolution expected from future radio telescopes and the next genera- 
tion of adaptive optics systems on large and extremely large telescopes. Further- 
more, anomalous flux ratios, astrometric perturbations, and time delay anoma - 



lies depend on different moments of the satellite mass function (|Keeton 2003 ) 



Therefore, a combination of techniques can help constrain both the slope and the 
normalization of the substructure mass function. 

6 Cosmography 

Cosmography is the measurement of the parameters that characterize the geom- 
etry, content, and ki nematics of the unive rse. Much progress has been achieved 
in recent years (e.g., iKomatsu et al. 20091 ). heralded as the era of precision cos- 



mology. However, some of the fundamental parameters need to be measured 
even more accurately if one wants to discriminate between competing theories. 
For example, the equation of state of dark energy w and its evolution with cos- 
mic time are essential ingredients to understand the nature of this mysterious 
phenomenon. 

Strong lensing is a powerful cosmographic probe, as it depends on cosmological 
parameters in two ways. Firstly, the time delay equation (and the lens equation) 
contain ratios of angular diameter distances. Therefore, within the context of a 
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model for the lensing potential, measurements of time delays or mass act as stan- 
dard rods, in a similar manner as the acoustic peaks of the power spectrum of the 
cosmic microwave background. Cosmography based on this concept is described 
in Sections 16. H 16.21 and 16. 3i Secondly, the optical depth for strong lensing de- 
pends on the number and redshift distribution of deflectors and therefore on the 
growth of structure and on the relation between redshift and comoving volume. 
Thus, given a model for the lensing cross section, and a model for the evolution 
of the population of deflectors, one can do cosmography from lens statistics. This 
approach is described in § 16.41 

6.1 Time delays 

Consider a galaxy lensing a time-variable source like a quasar or a supernova. 
Under the thin lens approximation, multiple images will be observed to vary 
with a delay which depends on the gravitational potential as well on a ratio 
of angular diameter distances (Equation [1]) . The ratio of angular diameter dis- 
tances is mostly sensitive to the Hubble Constant Hq (hereafter h in units of 
100 kms^^Mpc^^). However, time delays also contain non-negligible informa- 
tion about other cosmological parameters, especially if one considers a sample o f 



Coe fc Moustakas 20091 ) 



deflectors and sources spanning a range of redshifts (e.g. 
Therefore, although it is convenient to think in terms of the Hubble constant as 
the primary parameter, time-delays provide constraints in the multidimensional 
cosmological parameter space. When combined with other cosmology probes like 
the CMB power spectrum, time-delays are very effective at breaking degeneracies 
such as that between Hq and w (Figure [T3]) . 



Prom a practical point of view, cosmography with time-delays can be broken 
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into two separate problems: measuring time delays and modeling the lensing po- 
tential, including matter along the line of sight. Uncertainties in these two terms 
dominate the error budget and they are independent. Therefore, in order to mea- 
sure Hq to 1% accuracy from one lens system one needs to know both quantities 
with sub percent accuracy. Or, for a sample of N lenses, one needs unbiased 
measurements with approximately half ^/NYo uncertainty on both quantities. 
6.1.1 Measuring time delays Measuring time delays requires properly 
sampled light curves of duration significantly longer than the time-delay between 
multiple images. Once an approximate time-delay is known, the measurement 
can generally be refined by adapting the monitoring strategy, e.g. with dense 
sampling triggered after an event on the leading image. Typical time delays 
for galaxy lens systems are in the range weeks to months (with tails on both 
ends out to hours to years) and minimum detectable amplitudes from the ground 
are of order ~5%, limited by photometric accuracy for crowded sources and 
microlensing (see § 17. Sp . Thus, accurate time-delays typically require several 
seasons of dedicated monitoring effort. 



After the first "heroic" campaigns of the nineties and early 2000 (see lSchneider. Kochanek fc Wambsganss 2006 



for a review), which yielded of order 10 time-delays, several groups are now trying 
to take this effort to the next level with the help of queue mode scheduling and 
robotic t elescopes. A re cent summary of published time-delay measurements is 
given by iJackson (2007l ) . Tw o new time-delays have been published since then 



(J1206+4332 and J2033-4723; 



Paraficz. Hiorth k EHasdottir 2009 



Vuissoz et al. 20081 ). 



Taking the published time-delay uncertainties at face value, the present sample 
contributes to the error budget on Hq a little less than 1%. As I will discuss in 
§ El time-domain astronomy is a rapidly growing field and it is likely that many 
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of the logistical problems faced by time-delay hunters so far will be solved in the 

next decade. 

6.1.2 Determination OF THE LENSiNG POTENTIAL We now turn to errors 
associated with the local lensing potential, under the single screen approximation 
(matter along the line of sight and associated uncertainties will be described 
in ^ I6.1.3p . At fixed image configuration, time-delays depend to first order on the 
effective slop e of the mass distributio i i in the annular region between the multiple- 



images (see 



Saha fc Wilhams 20061 : 



Schneider. Kochanek Wambsganss 20061 . 



and references therein for discussion). For generic power-law models, at fixed 
lensing observables, the inferred Hq scales as Ho(7') ~ (7' — l)-ffo(7' = 2). For 
many systems, especially doubly imaged point sources, the lensing potential is 
highly uncertain and dominates the error budget. Unaccounted uncertainties in 
the mass model are the main culprits for the report ed discrepancies between t ime- 



Treu Koopmans 2002bl ^. 



delay determinations of Hq as large as ~30% (e.g. 

It is clear that some additional information is needed to bring the error budget 
on the lens modeling in line with that from time-delays. One approach con- 
sists of asserting some prior knowledge of the mass distribution in the deflectors 
and applying it to the analysis of a sample of systems. Since the effective slope 
is poorly constrained by lens data for point-like sources without additional in- 
formation, t he results depend critically on the prior. Following this approach, 



Oguri (20071 ) modeled 16 systems with power-law models assuming a Gaussian 
prior on 7' centered on 2 and width 0.15 obtaining h = 0.68 it 0.06 it 0.08 (the 
large systematic error attempts to reflect the large dispersion from system to sys- 
tem; however, it may also be due to the inclusion of systems with questionable 
redshift, time-delay or embedded in a complex cluster potential, which carries 
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substantial additional modelling uncertainties). The prior on 7 is plausible but 

not strictly justified, since there are no independent measurements for the sam- 
ple. For example, just changing the mean of the prior to (7') = 2.085lQQ^g it 0.1 
as found for the SLACS sample, would increase the estimate of Hq by 8%, with 
an additiona l systematic uncertainty of 10%. A very similar approach is that by 



Coles (20081 ). who imposes geometric priors to his pixelized mass reconstructions 
and obtains h = 0.71+^-^^ from 11 systems. Although it would be useful to draw 
samples from the lColes (20081 ) prior and measure the effective distribution of 7', 
it appears that his smoothness a nd steepness constraints create a distribution of 



effective slopes similar to that of jOguri (20071 ). explaining the agreement. These 



results are encouraging. However, they illustrate the challenge of reaching 1% ac- 
curacy using this methodology. One needs to have sufficient external knowledge 
of the distribution of mass in the sample of galaxies with measured time-delays 
to construct a sufficiently accurate prior. 

A more promising approach is to extract additional information for the very 
systems with measured time delays using ancillary data in addition to those 
available for the multiply-imaged point sources. In Bayesian terms, this means 
making the likelihood more constrainin g so as to reduce the relative impo rtance 



of the prior. Following this approach, 



Wucknitz. Biggs fc Browne (2004i ) mod- 



eled the extended radio structure around th e lensed quasar in B021 8+357 to 



Koopmans et al. (2003bl ) modeled 



infer 7' = 1.96 ± 0.02 and h = 0.78 ± 0.03. 
B 1608+656 using the measured stellar velocity dispersion and the HST images 
of the lensed host galaxy to measure 7' and infer h = 0.75^QQg ± 0.03, fixing 
= 0.3 and J^a = 0.7, and negelecting uncertainties due to the mass-sheet 
degeneracy, discussed in the next section. A recent analysis of improved Keck 
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and HST data of B1608+656 bv iSuvu et al. (201Q ) using more general pixellated 

models for the potential and the sou rce, infers h = 0.706ib0.031, for the same cos- 



mology as 



Koopmans et al. (2003b ). including uncertainties related to the mass- 



sheet degeneracy. This results shows that modeling errors can be reduced to a few 
percent per lens system, if sufficient observational constraints are available. If the 
other cosmological parameters are allowed to vary, one obtains the constraints 
shown in Figure [T3j The information fr om time delays is par ticularly powerful 
when combined to the WMAP5 results (jKomatsu et al. 20091 ). improving from 
h = 0.74^^^14 and w = -l-OdtlH to h = 0.697l|j:gf and w = -0.94l°;J^, for a 
flat cosmology. The results from a single lens are comparable with those from the 
local distance ladder r nethod (h = 0.742 zt 0.036 and w = —1.12 it 0.12 in combi- 



nation with WMAP5; iRiess et al. 20091 ) in terms of precision, although they are 
based on completely different physics and assumptions, and subject to different 
systematic errors. 

6.1.3 Mass along the line of sight and the mass-sheet degeneragy 
The flnal and perhaps limiting factor at this point is the uncertainty due to the 
unknown distribution of mass along the line of sight, i.e. deviations from the sin- 
gle screen approximation. On the one hand massive galaxies are typically found 
in groups. Group members and the common group halo contribute additional 
shear and convergence at the location of the main deflector. On the other hand, 
the "cone" between us the observer and source may be over or underdense, thus 
perturbing the time-delays with respect to those expected in a perfectly smooth 
and isotropic universe. Both effects can be thought to flrst order as equivalent to 
adding an external convergence Kext at the location of the deflector (which can 
be negative if the line of sight is underdense). 
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Due to the mass-sheet degeneracy, Kgxt is undetectable from dimensionless 

lensing observables. However, if we ignored its presence and make the standard 

assumption of vanishing convergence away from the lens to break the mass-sheet 

degeneracy, we would infer a biased value of Hq by a factor 1/(1 — Kext) (e.g. 



Schneider. Kochanek fc Wambsganss 20061 ) 



Independent measurements of mass, such as stellar velocity dispersion, help 
break the degeneracy because they constrain the local mass distribution. An un- 
known Kext leads to an overestimate of the lensing mass, and therefore alters the 
inferred 7' from comparison with kinematics, counterbalancin g the effects on Hp 



but not exactly. Mea, s urements of the local e nvironment (e.g, 



Fassnacht et al. 20061 . 



Auger et al. 2007 



Momcheva et al. 20061 ) also help, although the limiting fac- 



tor is the precision with which mass can be associated with visible tracers. A 
third approach consists of infer ring the distribution of effective Kcxt from high res- 



olution numerical simulations (jHilbert et al. 20071 ). The challenges of this third 



approach are producing realistic simulations at kpc scales relevant for strong lens- 
ing and understanding the selection function of the observed samples well enough 
to select simulated samples in the same way. In the case of B1608+656, the total 
uncertainty can b e brought to 5% using a combination of the three approaches 



(jSuvu et al. 2010l ). Analyzing a number of systems in similar detail will help 



uncover whether there are any residual significant biases. 
6.2 Lenses as standard masses 

Lensing studies indicate that the ratio fsiE between stellar velocity dispersion 
measured within a standard spectroscopic aperture and the normalization of 
the best fit SIE model ctsie is close to unity (1.019ib0.08 for the SLAGS sam- 
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pie for a concordance cosmology iBolton et al. 2008b ) , consistent with our gen- 
eral understanding of the mass distribution of early-type galaxies in the local 
Universe. If fsiE is known sufficiently well - independent of cosmology - lens 
galaxies could effectively be used as standard masses plugging measurements 
of Einstein radius and stellar ve locity dispersion into the SIE version of Eq. [9] 



(jGrillo. Lombardi &: Bertin 20081 ). Note that Hq cancels out in the ratio of an- 
gular diameter distances. Unfortunately, our current understanding of the mass- 
structure of deflectors and of the distrib ution of matter along the line of sig ht is 



not sufficient for accurate cosmography (jSchwab. Bolton &: Rappaport 2010l ). In 



some sense, the situation is similar to that of time-delay cosmography, and sim- 
ilar methodologies could be applied to overcome the limitations. The advantage 
of this method over time-delays is that it can be applied to any lens regardless 
of the presence of a variable source. The disadvantage is that the sensitivity of 
the angular diameter distance ratio on cosmological parameters is weak. 

6.3 Compound lenses 

The Einstein radius of a gravitational lens depends on the mass enclosed and on 
ratios of angular diameter distances. For systems with multiple sets of multiple 
images, such as SDSSJ0946-I-1006 (Figure [TU|) . one can solve for both the mass 
distribution and cosmography, provided that enough information is available to 
constrain the distribution of mass in the region between the Einstein rings. An 
additional complication is given by the mass associated with the inner ring, which 
acts as an extra deflector, making these syst ems compound lense s for the back- 



ground source responsible for the outer ring. 



Gavazzi et al. (2008i ) calculate that 



a sample of 50 systems like SDSSJ0946-I-1006 - expected for future large lens 
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surveys - should constrain the equation of state of dark energy w to about 10% 

precision. As in the cases discussed above, the issue is whether systematics asso- 
ciated with modeUng the deflector itself or the structure along the line of sight 
can be controlled with sufficient accuracy. High quality spatially resolved kine- 
matic information should help constrain the mass model of the main foreground 
deflector and of the inner ring. 

6.4 Lens statistics 

For a given source population, the fraction of strongly lensed systems (i.e. the 
optical depth) depends on the cross section of the deflectors and on the abun- 
dance of deflectors. Thus, measuring the abundance of strongly lensed systems 
constrains the intervening cosmic volume. This is the essence of lens statis- 
tics as a tool for cosmography, although quantities such as the distribution of 
Einstein radii and source redshifts also contain cosmographic information Note 
that lens-driven surveys are not nearly as sensitive as source-driven surveys (see 



Schneider. Kochanek &: Wambsganss 20061 . and references therein for a theoreti- 



cal description). 

The state of the art of thi s cosmographic application is the analysis of 11 



CLASS and 16 SQLS samples (|Chae 200?! . 



Qguri et al. 2 0081). w hich yield rather 



-0.3 



Qguri et al. 2008l l 



weak bounds on cosmological parameters (e.g., w = — 1.1±0.6_q 5 
Even though precision can certainly be improved by increasing sample size, the 
ultimate limit is set by systematic uncertainties. Accurate cosmography from 
strong lensing statistics requires accurate knowledge of: i) the mass-structure 
and shape of deflectors to compute cross-sections; ii) the contribution to the 
cross-section from large scale structures; iii) the number density of deflectors; iv) 
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the source luminosity function; v) the survey selection function. Quantities i to 

iv need to be known as a function of redshift. In conclusion, lens statistics poses 

three additional challenges (iii-v) over those in common with other cosmographic 

applications. 

6.5 Is lensing competitive? 

The studies mentioned in this section show that cosmography with strong-lensing 
gives results in agreement with independent probes, reinforcing the so-called con- 
cordance cosmology. However, the ultimate test for a method is when it breaks 
new ground in terms of precision, and the result is then confirmed independently. 
In my view, time-delays are the cosmographic application that stands the best 
chance of doing this for three reasons. First, two out of three major problems 
(time delay measurement and local mass model) have been solved and progress 
on the third (external convergence) is being made. Secondly, the inferred con- 
straints are well suited to break degeneracies inherent to other methods such as 
the CMB power spectrum. Thirdly, time-delays can be measured for a number 
of lenses using relatively small ground based telescopes or will come for free from 
future synoptic telescopes. Fourthly, the method is completely independent of 
the local distance ladder method and therefore provides a valuable independent 
test on its systematic uncertainties (like calibration and metallicity dependency 
of the cepheid-luminosity relation). Lenses as standard masses and compound 
lenses seem to be valuable cosmographic tools if they can be applied efficiently 
with limited observational resources, perhaps "piggy-backing" on other studies 

What is certainly very exciting and unique is the "inverse" application of cos- 
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mographic applications: learn about galaxy structure and evolution, on the basis 

of accurate cosmography from other probes. As mentioned in § HI time-delays, 

the combination of lensing and dynamics, and compound lenses have all been 

demonstrated to provide unique insights into the structure of distant galaxies, 

which cannot be obtained in any other way. This is also true for lens statistics, 

which can be used to determine the growth of the galaxy mass function in a 

unique way, once the m ass structure of each galaxy is un derstood from the other 



means discussed above f|Chae 200l iMitcheU et al. 20051 ') 



7 Lenses as Cosmic Telescopes 

In a typical galaxy-scale strong lens system, the background source is magnified 
by an order of magnitude. Exploiting this effect, lensed galaxies at intermediate 
and high redshift can be studied with the same level of detail as non-lensed galax- 
ies in the local universe (§ 17. ip . Furthermore, the host galaxies of bright active 
galactic nuclei are "stretched away" from the wings of the point spread function, 
enabling precise measurements of their luminosity and size, and ultimately of 
the cosmic evolution of the relation between host galaxy and central black hole 
(§ 17. 2p . Finally, microlensing by stars provides us with unique spatial informa- 
tion on the scale of the accretion disk, which is orders of magnitudes smaller than 
anything that can be resolved from the ground at any wavelength (§ 17. 3p . 

7.1 Small and faint galaxies 

The resolution of HST and the sensitivity of radio interferometers mean that 
we know very little about the distant (z >> 0.1) universe on scales below ~ 1 
kpc. Indeed, even in the nearby universe {z ~ 0.1), large ground based surveys 
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such as SDSS do not provide much sub-kpc scale mformation. Yet, we know 

from the local volume that small and faint galaxies are an essential ingredient 

of the universe, acting as building blocks of more massive systems. Only with 

the aid of gravitation al lensing we can reso lve sub-kpc scales and determine the 



morphology and size (jMarshall et al. 20071 ). and kinematics of small galaxies as 



well as trace the location o f star formation a nd the pattern of chemical abun 



dances ( Riechers et al. 2008 



Stark et al. 20081 ). Furthermore, flux magnification 



enables de tailed spectroscop ic studies that would be prohibitive in the absence 



of lensing (|Stark et al. 20081 ). These pilot studies show that intrinsic properties 
can be robustly recovered via lens modeling. The rapid increase in the number 
of known lenses should soon provide the large statistical samples needed for high 
impact studies. 



7.2 Host galaxies of lensed active nuclei 

In the local universe, massive galaxies are found to harbor central supermas- 
sive black holes. Remarkably, the mass of the black hole correlates with kpc- 
scale properties of the host bulge, su ch as velocity dispersion, luminosity and 



stellar mass (e.g. iGiiltekin et al. 20091 ). This family of correlations has been in- 
terpreted as evidence that black hole growth and energy feedback from active 
galactic nuclei play an important ro le in galaxy formation and evolution (e.g. 



Hopkins. Murrav &: Thompson 2003 ). However, the physics of the interaction as 
well as the relative timing of galaxy formation and black hole growth are poorly 
understood. Although the local relations are an important constraint, observ- 
ing their cosmic evolution is necessary to answer some fundamental questions. 
Are the local relations only the end-point of evolution, or are they established 
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early-on? Which comes first, the black hole or the host bulge? 

It is challenging to answer these questions observationally. Direct dynamical 
black hole mass measurements can only be done in the very local universe. At 
intermediate and high-z redshift, one needs to rely on indirect methods such as 
the empirically calibrated relation with continuum luminosity and line width ob- 
served for type-1 active galactic nuclei (AGN). However, th e presence of bright 



luminous point s ources hampers the study of the host galaxy (|Jahnke et al. 20091 . 



Treu et al. 200'il ). Strong lensing helps by stretching the host galaxy of distant 



lensed quasars primarily along the tangential direction (Figure [T^ . Of course, the 
quasar is also magnified, but one generally wins because the surface brightness of 
the point spread f unction falls off more rapidly than linearly. Using this method, 



Peng et al. (20061 ) showed that the bulges of host galaxies of distant quasars are 
more luminous than expected based on the local relation, consistent with a sce- 
nario where bulge formation predates black hole grow th, at least for so me objects. 
Similar results have been found for non-lensed AGN f^Treu et al. 2007i ). However, 
without the aid of lensing, studies have to be limited to lower redshifts and lower 
luminosity AGNs. 

7.3 Structure of active galactic nuclei 

Understanding the physics of accretion disks and the regions surrounding su- 
permassive black holes is essential to explain the AGN phenomenon with all its 
implications for galaxy formation and evolution. However, the scales involved are 
extremely small by astronomical standards (for a typical 10^ Mq black hole, the 
Schwarzschild radius is ~ 3 • 10^^ cm, the broad line region is ~ 10^''^^'^ cm), 
and therefore impossible to resolve with conventional imaging techniques. 
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Microlensing is perhaps the only tool capable of probing the small scales of the 

accretion disk. The Einstein radius of a star of mass (Fig. corresponds to 

approximately 4 • lO-'^^y^Ms/M0cm~ O.OIa/M JM^PC when projected at the red- 



Schneider. Kochanek fc Wambsganss 20061 ). 



shift of a typical lensed quasar (z^ = 0.5; Zg = 2[ 
The inner parts of the accretion disk will be smaller than this scale and there- 
fore subject to microlensing, while the broad line region and the outer dusty 
torus should be largely unaffected. The characteristic timescale for variation is 
given by the microlensing caustic crossing time, typically of order years, although 
i t can be shorter for special redshift comb inations such as that of Q2237+030 



([Schneider. Kochanek fc Wambsganss 20061 ) 



Based on this principle, one can infer the characteristic size of the accretion 
disk as a function of wavelength. Long-light curves - where the gravitational time 
delay betwe en multiple imag es can also be determined - provide the most strin- 



gent limits (jKochanek 20041 ). but interesting in formation can also be obtainec 



from single epoch 



data on a statistical basis (e.g. 



Bate. Webster fc Wvithe 2007 



Poolev et al. 200% . 

The inferred absolute size of the accretion disk can be known up to a factor of 
order unity, which depends on {Ms) and on the relative transverse speeds between 
the stars, the deflector, and the source. However, the slope of the relation between 
accretion disk temperature and size is independent of that factor and can thus 
be determined more precisely. Curre nt results indicate that the accretion disk is 
approximately the size expected for IShakura fc Sunvaev flOTa) disks, althoug h 



discrepancies of order a factor of a few have been reported ( 



Poolev et al. 200?! ). 



Assuming that the size sc ales as A^'^'^, t? is found to be in th e range 0.5 — 1, whereas 



7] = 0.75 is expected for a 



Shakura fc Sunvaev f 19731 ) disk (lEigenbrod et al. 20081 : 
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Poindexter. Morgan fc: Kochanek 20081 ). Long wavelength data imply the pres- 
ence of a second spectral component, consistent with the hypothesis o f a dusty 
torus of size much larger than the microlensing scale (jAgol et al. 20091 ). 

These first exciting results are just the beginning, because very few light curves 
obtained so far are long enough to harness the full power of microlensing. With 
the rapid development of time-domain astronomy predicted for the next decade, 
multiwavelength monitoring campaigns of several years for tens of objects should 
become feasible (§[9l). 

7.4 Cosmic telescopes and human telescopes 

I have described how strong lensing provides a unique opportunity to study 
sources that are too faint or too small to be studies otherwise, from quasar host 
galaxies to microarcsecond size accretion disks. 

Unfortunately, the use of galaxies (and clusters) as cosmic telescopes is often 
more contentious than it should be. One frequent critique is that source recon- 
struction is difficult and inherently uncertain. This is a false perception. The brief 
discussion in § [2] and the references listed therein provide ample documentation 
that lens modeling is now a mature field with very well understood uncertainties, 
capable of delivering results that are well reproduced by independent analyses. 
Lens modeling at cluster scales is more complex due to the larger dynamic range 
in the data and the more inhomogeneous mass distribution. Nevertheless, robust 
results can be obtained also for clusters, provided that enough information is 
available. 

Another frequent critique is that surveys using cosmic telescopes are inefficient 
compared to blank fields because of magnification bias. This is true for sources 
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with number counts in flux units (dN/dF) flatter than F~^. However, when 

probing the bright-end of the luminosity function of any population ~ where 

number density falls off exponentially - lensing is just unbeatable: the brighter of 

any class of distant astronomical objects will inevitably be gravitationally lensed. 

Cosmic telescopes and blank surveys are complementary to fully characterize a 

source population and its physical properties. 

8 Searches for Gravitational Lenses 

The strong lensing applications covered in this article span a broad range of 
astrophysical phenomena, observational, and theoretical challenges. However, 
they all share a common limitation: the relatively small number of systems to 
which they can be applied. Although there are 200 systems known, they are not 
all suitable for all applications. Studies must rely on at most a few tens of cases 
to infer results of general interest. 

Fortunately, a number of large surveys are expected to take place in the next 
decade, providing an ideal dataset to mine for rare objects such as strong lenses. 
The challenge will consist in developing fast and robust algorithms to find new 
lenses, and then in mustering the resources and the brain power needed to follow 
them up and study them (§ [9]). 

Before I summarize some of the searching techniques, it is useful to establish a 
discovery "etiquette": what are the necessary and sufficient elements to identify 
a strong gravitational lens? Here are two necessary criteria: i) multiple images 
clearly identified; ii) image configuration reproduced by a "simple" model. The 
first criterion seems to me unavoidable, although it has not always been applied in 
the past. The second criterion is more subjective, but can be made quantitative 
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in the following way. Given our knowledge of the surface brightness distribution 

of galaxies and of the gravitational potential, is it more likely that the observed 

configuration arises from some random configuration (e.g., HII regions distributed 

along a cross-pattern, or two quasars with similar colors on the opposite sides of a 

galaxy) , or from strong lensing of a more common surface brightness distribution? 

It seems to me these two criteria are also sufficient. Additional criteria such as 

images having identical colors or spectroscopic redshift of deflector and source 

are desirable, but impractical for future surveys that may have high resolution 

images in just a single band, or limited capabilities for spectroscopic follow-up. 

8.1 Imaging-based searches 

Imaging-based searches can be divided into catalog-based and pixel-based. Catalog 
based searches look for objects in a lensing-like configuration. They are most 
effective at d etecting sharp niultiply-imaged feat ures such as multiply-imaged 
quasars (e.g. Ilnada et al. 2008i . lOguri et al. 2008i ). but they can also be used 
for exten ded sources, prov i ded the image separa tion is large enough for de- 
blending fjAllam et al. 20071 . iBelokurov et al. 2007i ). Pixel-based searches start 
from a set of pixels, and look for lensing-like confi gurations. Lenses ar e iden- 



Cabanac et al. 2007 ) or by 



tified on the basis of characteristic geometries (e.g. 

actually modeling every system as a possible lens ([Marshall et al. 20091 ). The 
pixel-based method is slower and more computationally intensive than catalog- 
based searches, but in principle can be used to push the detection limit to 
smaller angular separations, beyond the level where source and deflector can 
be deblended by general-purpose cataloging softwares. Visual searches can be 
considered as pixel-based, with the human brain as lens-modeling tool (e.g. 
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Jackson 20081 : iNewton. Marshall &: Treu 2009i ). Algorithms need to be tweaked 

to reach an optimal balance between completeness (false negative) and purity 

(false positive) appropriate for each dataset and scientific goal. The best algo- 

rithm s can currently achieve 90% completeness and purity se arching through HST 



data ( Marshall et al. 2005 



Newton. Marshall fc Treu 20091 ). Although some hu 



man intervention is still necessary, this breakthrough makes it feasible to search 
through future surveys of 1000 deg^ or more. 

Time-domain surveys all ow for a different imag e-based strategy: looking for 



variable resolved sources (jKochanek et al. 2006al ). At high galactic latitude 



lensed quasars are more common than contaminants such as pairs of variable 
stars. Pairs of non-lensed quasars can be distinguished on the basis of their light 
curves and colors, while lensed supernovae are a welcome contaminant (see §ED- 
A first application of the method to the SDSS Supernovae survey data show that 
the only known compelling lens candidate is recovered as a close pair of variable 
sources. Out of over 20,000 so urces, only a hand ful of false positives are found. 



suggesting a "purity" of ~20% (jLacki et al. 20091 ). This is encouraging, although 



more tests on wider and deeper data are needed to further improve the method 
in view of upcoming surveys. 

8.2 Spectroscopy-based searches 

Spectroscopic searches rely on identifying composite spectra with features com- 
ing from multiple redshifts. Follow-up high resolution information is then needed 
to identify the subset of events with detectable multiple- images, and to obtain 
astrometry for lens modeling. A strong advantage of the method is that lenses 
come with redshifts by construction. After the early serendipitous discoveries 
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(|Huchra et al. 19851 ). the method starte d to bear large numbers of lenses only 



with the SPSS spectro scopic database (jAuger et al. 200S 



2006 



Bolton et al. 2008al . 



Willis et al. 20061 ). The recent searches highlight the quality of spectro- 
scopic data as the key element for success. High signal-to-noise ratios are needed 
to identify faint spectral features, close-to Poisson limited sky subtraction is 
needed to reduce false positives, spectral resolution better than 100 kms~^ is 
need to resolve line multiplets, and wide wavelength coverage increases the red- 



shift range for the search. It is a testament to t he high quality of t. 



he SDSS 



database that the confirmation rate is ~ 60 — 70 % ([Bolton et al. 2008al ) , after a 
very strict initial selection (approximately 1/1000 SDSS galaxies are selected as 
a candidate for follow-up by SLAGS). 



9 Future Outlook 



9.1 Thousands of gravitational lenses 

Most of the applications listed in the previous sections are listed by sample size. 
An increase by one of order of magnitude in sample size is needed to make 
progress. Fortunately, there is a realistic opportunity to make this happen in 
the next decade, considering the typical yields for strong lens systems searches. 
For optical and near i nfrared imaging searches, yields are ~ 10 deg-2 at HST-like 
depth and resolution fiMarshall. Blandford &: Sako 2005i ). ~1 deg^^ at the best 



ground based conditions 



Cabanac et al. 20071 ). At radio wavelengths and 0^.^25 



resolution expected for the Square Kilometer Array ^Koopmans et al. 2009al ) the 



yield is ~1 deg~^. For spectroscopic surveys, the yield is ~ 10~^/spectrum. Thus, 
a 1000 deg^ HST-quality cosmic shear survey, all sky ground based surveys in the 
optical or radio, and a 10^ galaxy redshift survey should all be capable of yielding 
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~ 10, 000 strong gravitational lens systems, although with different properties. 

High angular resolution surveys will be critical for applications such as the study 

of small mass deflectors and of the substructure mass function. Time-domain 

surveys will have a built-in advantage for, e.g., time-delays and microlensing. 

Spectroscopic surveys will be advantageous for those applications that require 

redshift and velocity dispersions, such as the study of luminous and dark matter 

in the deflector. Several thousand strong lens systems from each of these search 

techniques is an ambitious, yet feasible, goal for the next decade. 

These massive undertakings will require large number of people and resources. 

As in many other instances, it is likely that such projects will require the joint 

efforts of a number of communities interested in diverse scientific questions. The 

unique capabilities of strong lensing make it very worthwile to design future 

surveys keeping in mind its requirements. 

9.2 The problem of follow-up 

Let us assume that 10,000 strong lens candidates have been found. What follow- 
up will be needed to extract scientific information? Images with resolution of 
order O'.'l are often key to prove the lensing hypothesis, and to construct detailed 
lens models and study the properties of the host and the source. If the resolution 
of the finder survey is not adequate, follow-up will be required. Current follow- 
up imaging typically requires an orbit of HST. JWST should gain in speed for 
most applications and be revolutionary for long-wavelength studies, such as flux 
ratio anomalies. For a subset of objects with suitable colors and nearby stars, 
high resolution imaging could perhaps also be obtained in a comparable amount 
of time with a 8-lOm telescope equipped with laser guide star adaptive optics 
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(LGSAO). Extremely large 30-m class telescopes (ELTs) with LGSAO should 

be able to gain a substantial factor in speed and resolution. Radio follow-up 

of extended sources at high resolution with VLA requires of order 1 hour per 

lens. The Atacama Large Millimiter/Submillimiter Array (ALMA) should be an 

improvement both in speed and resolution. Even following up a thousand lenses 

will thus require thousands of hours of telescope time, maybe a few hundreds 

with JWST, ELTs and ALMA. This may be feasible, but not trivial, making high 

resolution imaging a likely bottle neck. Multiplexing is unlikely to be an option 

given the rarity of these objects on the sky, although multiplexing with different 

astronomical targets is certainly a desirable option. Even higher resolution images 

(O'.'Ol) are within reach with extreme adaptive optics on extremely large telescopes 

and will certainly be beneficial for pushing some of the lensing applications. For 

example, that kind of resolution could push the detection of DM substructure in 

distant galaxies in the 10^ Mq regime typical of the least massive luminous Milky 



Way satel 



ites currently 



mown, where the discrepancy with theory is currently 



strongest (jKravtsov 20101 ). 



Spectroscopic follow-up to gather redshifts is a problem of possibly even greater 
magnitude, considering that redshifts for many of the sources can not be measured 

Qfek et al. 2006 ). For the fainter 



even spending hours on the largest telescopes ((( 
sources, photometric redshifts may be the only option. Coordination with redshift 
surveys - such as those proposed to measure baryonic acoustic oscillations - will 
help measuring redshifts as well as in spectroscopic searches, although they will 
also require high angular resolution follow-up. 

Monitoring campaigns of thousands of lensed AGNs are out of the question at 
the moment, but could be a natural byproduct of future synoptic surveys. Some 



58 Tommaso Treu 

of the most demanding time domain applications, such as detection of time-delay 

anomalies could be beyond the reach of ground bas ed monitoring tools and require 



a dedicated space mission ([Moustakas et al. 20081 ) . 



In parallel with discovery efforts, careful thought must be put into planning 
follow-up efforts. First, ways to extract as much information as possible from the 
discovery images themselves must be found. Second, follow-up efforts should be 
coordinated as much as possible with those of other science cases to find common 
paths and synergies. Last but not least, brain power could be another serious 
limitation. Currently, accurate and reliable lens models require several days of 
expert human brain activity. This will not be possible when samples will consist 
of tens of thousands of systems. 

9.3 Unusual lensing applications in an era of abundance 

I conclude with four examples of strong lensing applications that require very 
rare conditions and therefore need the large samples expected in the next decade 
to become viable. 

Lensed supernovae la are extremely valuable, because their standard luminos- 
ity constrains the absolute magnification and therefore breaks the mass-sheet 
degeneracy. For typical rates, we expect of order one could be found monitoring 
known lenses. However, a ground based time-domain survey covering most of the 
sky is expected to find of order a hundred lensed type las (Oguri &: Marshall, 
2009, in preparation). 

Compound lenses are potentially powerful cosmographic probes, but there is 



currently only one such system known at galaxy scales (jCavazzi et al. 20081 ) 



Thousand square degree field surveys at HST-like resolution should be able to 
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find ten s of systems like Sp SSJ0946+1006, potentially constraining w to the 10 



% level (|Gavazzi et al. 20081 ). 



Strong lensing is one of the few tools capable of measuring the mass of qui- 



escent black holes at cosmo^ 
properties of central images 



ogical distances, through thei r gravity affects the 



Mao. Witt k Kooomans 2001 



Detecting the cen- 



tral image - which is generally highly demagnified - is usually beyond reach with 



current instrumentation (see however 



Winn. Rusin &: Kochanek 20041 ). However 



this application may become practical with future facilities, especially at radio 
wavelengths where the contrast between deflector and source is more favorable. 
Finally, with futur e samples of 10^ lenses, rare exam ples of "catastrophes" 



should be identifiable (jOrban de Xivrv &: Marshall 20091 ). These are very special 



lensing configurations characterized by specific constraints on the gravitational 
potential and its deri vatives, and they only occur only for very specific source posi- 



tion and redshift (see 



Fetters. Levine fc Wambsganss 2001 



Schneider. Ehlers k Falco 19921 . 



for details). The identification of examples of catastrophes is interesting for two 
reasons. Firstly, catastrophes often lead to extreme magnification factors, up 
to ~ 100, making them extraordinary cosmic telescopes. Secondly, the unusual 
geometry of multiple images can give remarkably strong constr aints on the mass 



distribution of the deflector (jOrban de Xivrv Marshall 20091 ). 
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Acronyms 

1. CDM: Cold Dark Matter 

2. CMB: Cosmic Microwave Background 

3. DM: Dark Matter 

4. HST: Hubble Space Telescope 

5. IMF: Initial Mass Function 

6. JWST: James Webb Space Telescope 

7. SDSS: Sloan Digital Sky Survey 

8. SIE: Singular Isothermal Ellipsoid 

9. SIS: Singular Isothermal Sphere 

10. SLACS: Sloan Lens Advanced Camera (for Surveys) Survey 

Definitions 

1. Convergence: dimensionless projected surface mass density in units of the 
critical density. 

2. Deflector: the foreground galaxy responsible for the lensing potential. 
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3. Einstein radius: characteristic scale of strong lensing. For a circular deflec- 
tor it corresponds to the radius within which (k) = 1. 

4. Image plane: twodimensional map of the source emission as it appears to 
the observer after propagation through the lensing potential. 

5. Macrolensing: strong lensing producing image separations of order arcsec- 
onds, the typical scale of massive galaxies. 

6. Microlensing: strong lensing producing image separation of order of micro- 
arcseconds, the typical scale of individual stars. 

7. Millilensing: strong lensing producing image separation of order of milli- 
arcseconds, the typical scale of small satellite galaxies. 

8. Shear: dimensionless quantity that describes the local distortion of lensed 
images. 

9. Source: the background astronomical object whose light is being lensed. 

10. Source plane: twodimensional map of the source emission as it would appear 
to the observer in the absence of a deflector. 

11. Strong lensing: deflection of light from a background source by a foreground 
deflector strong enough to produce multiple images. 

Summciry Points 

1. Massive early- type galaxies are surrounded by dark matter halos spatially 
more extended than the luminous component. The fraction of mass in the 
form of dark matter inside the effective radius increases with galaxy stellar 
mass. 
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2. The total mass density profile of massive early- type galaxies is approxi- 
mately isothermal in the innermost ~10 kpc, i.e. the logarithmic slope 7' 
equals two within 10%. 

3. Precise gravitational time delays for a single system can be used to measure 
the Hubble Constant to 5% precision, provided that enough information is 
available to constrain the local gravitational potential and to break the 
mass-sheet degeneracy. Time delays break the degeneracy between h and 
w in the analysis of CMB data. Combining the constraints from the lens 
system B1608-F656 and those from WMAP5 yields h = Omitoml and 
w = —O.M^q'II assuming flatness. 

4. The host galaxies of distant luminous quasars appear to be underluminous 
in comparison with local galaxies hosting black holes of the same mass. 
This may indicate that in this mass range black holes complete their growth 
before their host galaxy. 

5. Microlensing results indicate that the size of accretion disks and its de- 
pendency on temperature is in broad agreement with the predictions of 



Shakura Sunvaev (19731 ) models. Moreover, mid-infrared microlensing 
studies are consistent with a presence of an unresolved dusty region, larger 
than the accretion disk. 



Future Issues 

1. How do luminous and dark matter density profiles depend on galaxy mass, 
type, and cosmic time? 

2. Are dark matter density profiles universal, as predicted by CDM numerical 
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simulations? 

3. Is the mass function of substructure in agreement with the predictions of 
CDM numerical simulations? 

4. How are density profiles and the substructure mass function influenced by 
the presence of baryons? 

5. Is dark energy the cosmological constant {w = —1)? If not, how does the 
equation of state evolve with cosmic time? 

6. How can we find and exploit larger samples of strong gravitational lens 
systems? 
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Figure 1: Sketch of the gravitational lensing geometry, courtesy of B. Brewer. 
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Figure 2: Examples of the most common configurations of galaxy-scale gravi- 
tational lens systems. A background source (top left) can produce four visible 
images (a "quad"; top right), an (incomplete) Einstein ring (bottom left), or 
two visible images (a "double"; bottom right), depending on the ellipticity of 
the projected mass distribution of the d eflector and on the re lative alignment 
between source and deflector (data from iMoustakas et al. 20071 . Image courtesy 
of P. Marshall). 
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Figure 3: Optical analogy to illustrate the gravitational lensing phenomenon. 
The optical properties of the stem of a wineglass are similar to those of a typical 
galaxy scale lens. Viewed through a wineglass, a background compact source 
such as distant candle (top left), can reproduce the quad (top right), Einstein ring 
(bottom left), and double (bottom rights) configurations observed in gravitational 
lensing and shown in Figure [2l Image courtesy of P. Marshall. 
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Figure 4: Einstein radius of a massive elliptical galaxy (top), a dwarf satellite 
(middle) and a star (bottom) as a function of deflector redshift for three choices 
of source redshifts (zg = 1, 2, 8). SIS models with velocity dispersion (T=300 and 
10 kms~^ are assumed for the elliptical and dwarf galaxies respectively. A point 
mass of one solar mass is adopted for the star. 
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Figure 5: Example of a gravitational lens model, from iBolton et al. (2008al . re- 
produced by permission of the AAS). The two left panels show the data before 
and after subtraction of the light from the lens galaxy. The smaller panels on 
the right show the predicted image intensity of the best fit lens model, residuals, 
and source plane reconstruction, for an SIE mass model (top panels) and a mass 
traces light model (bottom panels). In the panel representing the image plane 
(labelled SIE) the white line shows the critical line. In the panel representing the 
source plane (magnified by a factor of 2) the white lines show the inner and outer 
caustics. Note that the peak of the surface brightness distribution is located out- 
side the inner caustic and is therefore imaged twice, while the outer regions of 
the lensed sources go through the central region and therefore form an Einstein 
ring in the image plane. 
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Figure 6: Microlensing observed in the quadruply-imaged quasar PG1115+080 
{zs = 1.72). The lens galaxy (zd = 0.31) has been removed for clarity. Each 
panel is 4 arcseconds on a side. The bottom right panel (labelled Model) shows 
the expected image predicted from an SIS model of the deflector and an external 
shear term to account for the effects of a nearby group. The flux of image A2 
incre ased by over a fact or of four between June 2000 and January 2008 (Figure 



from 



Poolev et al. 2003 . reproduced by permission of the AAS). 
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Figure 7: Distribution of deflector (top) and source (bottom) spectroscopic red- 
shifts for the galaxy-scale gravitational lens systems discovered by the CLASS, 
COSMOS, SLACS, and SQLS surveys (see text for details). 
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Figure 8: Dark matter fraction inside the cylinder of projected radius equal 
to the Einstein radius as inferred from stellar population synthesis modeling of 
multicolor data and strong g ravitational lensing analysis of the SLACS sample 



(data from lAuger et al. 20091 ) 
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Figure 9: Mass density profiles o f lens galaxies inferred fro m a strong lensing and 



dynamical analysis (Figure from 



Treu &: Koopmans 2004 . reproduced by permis 



sion of the AAS). In addition to the mass associated with the stars (red line), 
the data require a more extended mass component, identified as the dark matter 
halo (blue line). Although neither component is a simple power-law, the total 
mas profile is close to isothermal, i.e. 7' = 2. The vertical dashed line identifies 
the location of the Einstein radius. 
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Figure 10: Double Einstein ring compound lens SDSSJ0946+1006. Left: color 
composite HST image (Courtesy of M. W. Auger). Note the foreground main 
deflector in the center, the bright ring formed by the images of the intermediate 
galaxy, and the fainter ring formed by the images of the background galaxy lensed 
by the two intervening objects. Right: Enclosed mass profile as inferred from 
the Einstein radii of the two rings (red solid points - the error bars are smaller 
than the points). The enclosed mass increases more steeply with radius than 
the enclosed light (solid blue line; rescaled by the best fit stellar mass-to-light 
ratio), indicating the presence of a more extended dark matter component. Even 
a "maximum bulge" solution (dotted blue line) cannot account for the mass at 
the outer Einstein radius. 
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Figure 11: Example of edge-on spiral lens system (z^ = 0.063) discovered by 
the SWELLS Survey. The multiply-imaged source {zg = 0.637) is visible in the 
optical HST discovery image and readily apparent in the Keck near infrared image 
where the effects of dust are minimized. The combined information at multiple 
wavelengths allows one to correct for dust and infer the stellar mass of the disk 
(Lnage credits: A. Dutton, P. Marshall, T.Treu). 
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Figure 12: Near infrared (2.2/xm) image of the gravitational lens system 
B2045 +265 taken with the adaptive optics system at the 10m Keck-II Telescope 
(from iMcKean et al. 2007i ) . For this kind of configuration, the flux of image B 



is expected to be equa. 



sence of substructure (IBradac et al. 20021 : 



to the sum of the fluxes of images A and C, in t he ab- 



Keeton. Gaudi fc Fetters 200.'^ ). The 



anomaly was originally discovered on the basis of radio images, ruling out mi- 
crolensing or differential interstellar medium scattering as alternative interpre- 

n 1 1 ' 

tation (jFassnacht et al. 19991 . IKoopmans et al. 2003a ). A satellite galaxy (G2) 
of the main deflector (Gl) is detected in this deep and high resolution image. 
A small mass located at the position of the satellite can explain the observed 
anomaly. 
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Figure 13: Illustration of cosmography with gravitational time delays. The panels 
show two and one dimensional posterior probability distribution functions for Hq, 
w, and assuming flatness. Red lines indicate limits from cosmic microwave 
background, blue lines represent limits obtained from a single gravitational lens 
system with measured time delays (B1608+656), while black lines represent the 
joint constraints. Note how the constraints from time-delays are almost vertical 
in Hq and therefore help break the degeneracy between w and Hq in the CMB 
data. Lensing constraints in the w — are broad and therefore not shown for 
clarity (Figure courtesy of S.Suyu; Figure from Suyu et al. 2010) . 
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Figure 14: Illustration of gravitational lenses as cosmic telescopes. Two-image 
lens system (HE 1104-1805) of a Zg = 2.32 quasar produced hy a = 0.73 
foreground galaxy. Panel a) shows the original data, panel b) shows the lensed 
host galaxy found after subtracting the deflector and quasar components of the 
best-fitting photometric model, panel c) shows the residuals from that photomet- 
ric model, and panel d) shows what the unlensed host galaxy would look like 
in a similar exposure after perfectly subtracting the flux from the quasar. The 
curves shown superposed on th e model of the host galaxy are the lensing caustics. 



(Figure from 



Peng et al. 20061 . reproduced by permission of the AAS) 



